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We investigate the application of a perturbative renormalization group (RG) method to random 
antiferromagnetic Heisenberg chains with arbitrary spin size. At zero temperature we observe that 
initial arbitrary probability distributions develop a singularity at J = 0, for all values of spin S. 
When the RG method is extended to finite temperatures, without any additional assumptions, we 
find anomalous results for S > 1. These results lead us to conclude that the perturbative scheme 
is not adequate to study random chains with S > 1. Therefore a random singlet phase in its more 
restrictive definition is only assured for spin- 1/2 chains. 



Low-dimensional random quantum magnetic systems have become an object of increasing interest in recent years. 
Despite the apparent simplicity, quantum spin chains show a wealth of physical properties which have attracted the 
attention of both theoretical and experimental physicists. It is generally accepted that quantum antiferromagnetic 
spin-1/2 chains in the presence of disorder exhibit, at low temperatures, a random singlet (RS) phase characterized by 
power law behavior of thermodynamic quantities [|J0]- This has, in fact, been observed in several experiments 
Theoretically, the new phase was revealed by a real space renormalization group (MDH) approach jj] . These results 
have motivated many authors to study the physical properties of other random antiferromagnetic quantum (RAQ) 
chains with S > 1. In this case there are additional new problems, as for example, what happens to the Haldane 
gap of integer spins chains || in the presence of disorder. In fact, the possibility of the energy gap be suppressed by 
disorder, driving the system to a RS phase has become a very attractive matter for investigation. In particular, the 
5 = 1 RAQ chains has been exhaustively studied by several methods |)] |l6|. Many authors have employed extended 
versions of the MDH scheme to explore the ground state properties of the model. In the course of these investigations, 
some authors suggested the existence of random singlet phase in the spin-1 chain in the strong disorder regim e p ft3[ . 
On the other hand, studies based on exact diagonalizatio n Jl4| , quantum Monte Carlo (MC) simulation ]15| and 
density matrix renormalization group (DMRG) techniques ]16[ found that the Haldane phase is quite robust against 
randomness. 

In view of these controversial results and the fact that most calculations for spin-1 chains have been done at T = 0, 
we were motivated to extend the MDH analysis to finite temperature. Our goal was to generalize the MDH recursion 
relations to study the thermodynamic behavior of chains with arbitrary spin size. We performed careful numerical 
procedures to iterate our generalized MDH equations at finite temperature. The anomalous results we have obtained 
for the free energy and are presented in this communication indicate that the MDH approach breaks down for chains 
with S > 1. Therefore, we argue that it is impossible to conclude on the existence or not of the RS phase in the S > 1 
RAQ chains, as proposed in previous zero temperature studies using the MDH method p|-|l3|]. Actually, within the 
MDH approximation the existence of the RS phase is assured only for the spin-1/2 chains. 

The random magnetic systems are described by the Heisenberg Hamiltonian: 

L-i 

H = 2_j JiSi.Si+i, (1) 

where Ji > is the nearest neighbor antiferromagnetic interaction and Si are quantum spin operators. The exchange 
couplings are random variables distributed according to a given probability distribution Pj(Jj,f2), with a cutoff Q. 
The MDH method consists in eliminating the pair of spins with the strongest coupling ( J 2 = £1) in the random chain 
by considering the interaction ( J\ and J3) with the neighboring spins of this pair as a perturbation (see Fig. Q). 
The Hamiltonian for the strongly coupled pair of spins S 2 and S3 is given by 

H = J 2 S 2 .S 3 . (2) 
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These spins are weakly coupled to the neighbors via 



Hi — J\S\.S2 + J3S3.S4. (3) 

We calculate the free energy of the system formed by S%, S2, S3 and 64 through perturbation theory in second order 
of Hi (see ref. After some simple but extensive calculations we come with the results: 

F' = F - ±[S(S + l)] 2 (Ji 2 + 4)V s m) (4) 

with 

25 

F = -S(S + 1)^-4 ln£(2i + 1) exp[--i(i + (5) 



and 



J' = -S-(S + l)^W s (/30), (6) 



where /3 = 1/kgT and the functions V and W are given by 



v n _ (gg+ !) 2 - £-fo(2* + l)ej^[l - W + l)y] 
W 4S(S+l)£ 2 f (2z + l) e -^+ 1 )*> 



and 



(25 + I ) 2 - E'fo(2» + l)e-^ i+1 )»[l + + l)y] 
4S(S + 1) £ 2 £ (2i + l)e-**(*+*)v 



w.(y) = - \^:^s Z. . - _^:r — ■ w 



with 



Equations (U)-® are the generalization of the MDH equations for the thermodynamic properties of RAQ chains with 
arbitrary spin size. For S=l/2 we recover the original MDH recursion relations [Q. The S—l case has already been 
studied in the limit of very low temperatures flo[| . 

When the decimation process is carried out for a spin-1/2 chain the distribution Pj(Ji,£l), independently of the 
form of the original distribution of antiferromagnetic bonds, rather quickly approaches a fixed form singular at J = 

§, 

"—in ™ 

The fixed form of Pj(Ji, fi) determines the low-temperature behavior of the thermodynamic quantities. The expo- 
nent a gives rise to logarithmic corrections in the thermodynamic functions. 

The great merit and the reason for the success of the MDH theory H , was to show that the characteristic RS phase 
behavior of spin-1/2 chain is universal. More specifically, it is independent of the original distribution of exchange 
couplings. The MDH approach represented a considerable progress with respect to previous methods |5|||[l7j which 
required an initial distribution of bonds with a particular form. Actually, in order to reproduce the power law 
dependence of the low temperature thermodynamic quantities observed experimentally, the previous models required 
a initial distribution already with a power law behavior P,|l7|] . 

We turn now to the numerical simulation of chains with spins S=l and S=3/2. Each chain is composed of N 
spin-S objects with periodic boundary conditions (with N= 50,000). We carried out the averages over 10 different 
configurations for each Pj(Ji, f2) to obtain the free energy of the system. We start by choosing the exchange coupling 
from an uniform distribution Pj(Ji, f2) given by: 

f t^a for A < Ji < 1, 

Pj(j i) n) = i !- A - ' (11) 

Jy ' ' I otherwise V ' 
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The parameter A represents the strength of the initial disorder of the couplings. The case A = represents the strong 
disorder limit because the distribution becomes extremely broad in a logarithmic scale; namely, infinitesimally weak 
bonds appear. The weak disorder regime is represented by finite values of A. The distribution, in this case, presents 
a gap A. 

We have studied the flow of the initial coupling distribution, Eq. (jll]), for the spin-1 chain for the following cases: 
A = 0, A = 0.05, A = 0.1, and A = 0.2. We found that for any A, successive elimination transformations give rise to 
weaker and weaker couplings as the cutoff 51 decreases, independently of the temperature. For sufficiently small 51, i.e., 
after a sufficient number of eliminations, the distribution Pj(Ji, 51) becomes peaked at J,; = and, as in the T=0 case 
- |l3[ , can be approximated by a power-law. We should point out that the convergence to a power-law depends on 
the value of A. It is faster in the case of the gapless (A = 0) initial distribution. As in the case of spin- 1/2 chains § 
the power-law exponent obtained for the S=l chain is non- universal: it depends on temperature and cutoff, although 
weakly. In Fig. we present, for completeness, the power-law behavior of Pj(Ji, 51) with exponent a = 0.38 obtained 
from an initial distribution with A = 0, at a fixed temperature ksT = 0.1051. 

Now we focus on the free energy per spin and the specific heat of the spin-1 chain. In Fig. || we show the plot of the 
specific heat versus temperature (given in units of 51). Despite the power- law dependence of Pj(Ji, 51) at low energies 
described above, we find a non physical behavior in the free energy of the system. For intermediate temperatures, 
the free energy gives rise to a negative specific heat (see Fig. ||). Nonetheless, at very low temperatures the free 
energy and thus the specific heat is well-behaved, being described by a power-law (see inset of Fig. ||). These very low 
temperature features are similar to the random singlet phase of the disordered spin-1/2 chain ffl. Moreover, at higher 
temperatures, the free energy curve is again well-behaved and gives rise to a physically meaningful specific heat. 

The presence of the factor 2S(S + l)/3 in Eq. (||), which defines the renormalized coupling J', plays an important 
role in our discussion. In the case of the spin-1 chain it is 4/3. At low temperatures (fceT < 0.3751), the function 
Wi(pfl) is larger than 3/4, consequently, when the couplings neighboring the cutoff 51 are not small enough, couplings 
larger than those eliminated may be generated (see Eqs. (||) and (||)). In this low temperature regime, although the 
cutoff decreases rapidly, large couplings are generated with non negligible probability. When the couplings between 
the strongly coupled pair of spins and its neighbors are not sufficiently weak, the perturbative scheme fails, giving 
rise to the spurious result observed in the free energy. At higher temperatures (fc^T > 0.370) the MDH method gives 
correct results for the thermodynamic properties pp| . There the function Wi(/3Q) < 3/4 cancel the coefficient 4/3 of 
J{ (see Eq. (||)), thus validating the present approach, independently of the starting distribution. In addition, when 
one considers distributions with a finite gap (weakly disordered distributions), the breakdown of the MDH formalism 
is most pronounced. The probability of obtaining, at low temperatures, the product J\ J3 < 3/451 2 depends on the 
value of the gap. By increasing A, the initial distribution will not generate infinitesimally weak couplings to ensure 
that the product is sufficiently small. 

In the spin-3/2 chain we performed the decimation process by starting with a uniform gapless probability distribu- 
tion. Despite the larger coefficient in Eq. (||), Pj(Ji,Q) still develops a power-law behavior as the cutoff 51 becomes 
sufficiently small. However, this convergence is much slower than that for a spin-1 chain. Now the couplings J' which 
are larger than the decimated out coupling 51 are statistically predominant. The effect of the breakdown of the MDH 
method in the free energy of the S — 3/2 chain is more dramatic yielding a negative specific heat even in the very 
low temperature limit. 

These results has led us to try to improve the MDH transformation in two directions: first we considered the 
decimation of larger spin clusters. In order to gain an insight on whether or not this procedure would work, it was 
sufficient to carry out the decimation at zero temperature. The intent was to investigate whether the perturbation 
scheme, using larger clusters, would decrease the factor appearing in the expression for the new coupling J'. We 
have carried out numerically the decimation process for RAQ chains with both spin S — 1/2 and S = 1, considering 
clusters with 2, 4, 6 and 8 spins. In our numerical calculations we considered, for simplicity, Jj = 51 = 1, where the 
J[s stand for the bonds between neighboring spins in the cluster. The results for the case of a cluster with only two 
spins can be directly compared with the value 1/2 and 4/3 of the original MDH procedure (see Eq. (||)) in the S=l/2 
and S=l RAQ chains, respectively. We list the results in Table ffl. In the case of a spin-1/2 chain the results indicate 
that the MDH scheme is always improved. However, the resulting values for larger clusters in the spin-1 chain makes 
clear that the problem remains unsolved. 

We also tried to improve the MDH scheme by considering terms of higher order in the perturbation expansion, at 
zero temperature. In that case, our purpose was to investigate whether the higher order corrections would decrease 
the factor 2S(S + l)/3 in Eq. (||). We have thus calculated the T — renormalized coupling, for the spin-1 chain, up 
to third order in 51, which is given by 
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J, = |] JiJ3 + i (Ji2j3 + JiJ 3 ] (12) 

Once again, we come upon a result which reinforces that the MDH formalism is not suitable to describe the low 
temperature properties of the RAQ chains with S > 1 (defined in Eq (0)). Moreover, in the fourth-order correction to 
the renormalized bond biquadratic exchange couplings are generated. In this case we should consider these interactions 
already from the beginning in the Hamiltonian model. Finally we call attention that starting with a power law 
distribution, the MDH procedure works very well. In this case, only pairs of spins weakly coupled are present in the 
chain and the probability of the generated coupling being larger than the one eliminated is statistically negligible. 
The singular form of Pj(Ji, f2) dominates the thermodynamic properties which are also described by power-laws. Of 
course this is in strong contrast with the robust and universal character of the RS phase for the spin-1/2 chain, which 
is indeed a new fixed point 0,^1- 

In summary, at zero temperature there is not a safe criterion to discuss the possibility of the occurrence of the RS 
phase in the RAQ chains. In fact, it was necessary to extend the MDH scheme, for S > 1, to finite temperature to 
establish the failure of the perturbative approach. In view of the results we have obtained, we strongly claim that the 
MDH procedure is not adequate to study the RAQ chains with S > 1 defined in Eq. (Q). Moreover, the existence of 
a RS phase at low temperature is a property assured only for spin-1/2 RAQ chains. 
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TABLE I 



cluster size 


spin- 1/2 


spin-1 


2 


0.500000 


1.333333 


4 


0.377992 


1.881373 


6 


0.312648 


2.597776 


8 


0.270240 





TABLE I. Renormalized coupling for chains with 5 = 1/2 and S = 1 obtained with clusters of size 2,4,6 and 8. We have 
taken Ji — Q = 1. 



FIGURES CAPTION 



FIG. 1. Spins and coupling constants involved in the elimination transformation (a) to (b). 

FIG. 2. Power-law behavior of Pj(Ji, fi) in the low-energy limit, at UbT — O.lOfl Inset: the Log-Log plot of Pj(Ji, fi) yields 
the exponent a = 0.38. 

FIG. 3. Low temperature behavior of the specific heat of the random spin-1 chain for A = 0. Inset: the power-law behavior 
with exponent equal to 0.77 
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